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Abstract 

We study the effect of inhomogeneities in the magnetic field on the thermal 
entanglement of a two spin system. We show that in the ferromagnetic case a very 
small inhomogeneity is capable to produce large values of thermal entanglement. 
This shows that the absence of entanglement in the ferromagnetic Heisenberg 
system is highly unstable against inhomogeneoity of magnetic fields which is 
inevitably present in any solid state realization of qubits. 
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1 Introduction 

1.1 Motivation 

It is well known that quantum entanglement [21 El plays a fundamental role in al- 
most all efficient protocols of quantum computation (QC) and quantum information 
processing 01 

Without entanglement which is the essential quantum ingredient of QC, any quan- 
tum algorithm which only uses the other property of quantum mechanics, namely 
the superposition property, can also be implemented on any physical system which 
allows superposition of states, i.e. classical linear optical devices. In any proposal for 
physical implementation of qubits, it is therefore of utmost importance to investigate 
the entanglement properties of pairs and collection of such qubits. Among the many 
proposals for physical implementation of qubits, those based on solid state devices 
seem to be promising as far as the crucial scalablity property is concerned. 
In one such proposal [S] a well localized nuclear spin coupled with an electron of a 
donor atom in silicon plays the role of a qubit which can be individually initialized, ma- 
nipulated and read out by extremely sensitive devices. In another proposal [7|IH1 1^1 [TTI]. 
the spin of an electron in a quantum dot plays the role of a qubit. Long decoherence 
time and scalability to more than 100 qubits are two of the important virtues of this 
scheme. 

In both schemes the effective interaction between the two qubits is governed by 
an isotropic Heisenberg Hamiltonian with Zeeman coupling of the individual spins, 
namely 



Actually the isptropic interaction is an approximation, since spin orbit coupling intro- 
duce perturbations which break this isotropy . A more complete hamiltonian would 



of the order of 10~^. Note that the only coupling in the interaction between spins that 
is controllable is J ^21; and the individual couplings between different components of 
spins denoted usually by J^, Jy and Jz can not be controlled separately and thus one 
can not adjust these parameters arbitrarily to enhance the entanglement in a given 
situation. 

This means that although studies of entanglement for different types of anisotropic 
interactions are very interesting theoretically (specially when infinite spin systems are 
treated which is the only case which yields valid results with regard to quantum phase 
transitions pSl); they may not be of much practical relevance to concrete physical 
realization of qubits. 

In this paper we ignore the anisotropic perturbations due to both their smallness and 
due to the fact that strategies have been invented to cancel such anisotropics |14j . 
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Due to their smallness, they may only introduce minor changes in any result derived 
for the isotropic case. 

On the other hand in any solid state construction of qubits, there is always the 
possibility of inhomogenous Zeeman coupling ^1 . Solid state heterostructures 
are usually inhomogeneous and magnetic imperfections or impurities are likely to be 
present leading to stray magnetic fields. Indeed it is one of the main challenges in 
this proposal to construct identical qubits J7j . Constructing nearly identical devices 
in semiconductor technology has always been difficult and is still difficult, e.g. a very 
small temperature or strain difference in the substrate produces differences which al- 
though may not be significant for the classical semiconductor technology will certainly 
be important for the quantum technology ^Zj. Besides these unwanted effects, there 
are schemes like parallel pulsed schemes ^Sl) in which both a localized and hence 
inhomogenous Zeeman coupling and exchange interactions are employed to expedite 
manipulation of qubits. 

In view of the above it is desirable to consider a two qubit system in an inhomogenous 
magnetic field and study the entanglement properties of this system in detail. 
At extremely low temperatures such a qubit system may be assumed to be in its 
ground state. Thus it will be desirable to study the entanglement properties of the 
ground state. 

However a real physical system is always at a finite temperature and hence in a mix- 
ture of disentangled and entangled states depending on the temperature. Therefore 
one is naturally led to consider the thermal entanglement of such physical systems. 
In summary we mean that the thermal entanglement of finite systems has more rele- 
vance to the problem of initialization of quantum computers jl9j than to the problem 
of quantum phase transitions which requires a study of infinite size systems. 

1.2 A brief account of previous works 

Thermal entanglement in a two qubit Heisenberg magnet with the Hamiltonian 



was first studied by Nielsen [21] who showed that in the ferromagnetic case (J < 
0) no entanglement exists but in the antiferromagnetic case (J > 0) entanglement 
appears below a threshold temperature Tc. Since then many other systems have been 
investigated. 

There is now a vast literature on this subject and for clarity it is better to separate 
them into two categories, namely those ^22 , 23 , 24tl25|[26] which study by analytical or 
numerical methods infinite spin chains with at times particular attention to quantum 
phase transitions and those which study few, mostly two, spin systems. In our opinion 
one can not draw valid results for quantum phase transitions by studying a two spin 
system, and these types of studies are useful in other contexts, e.g. the problem of 
initialization of a quantum computer as described above, provided they start with a 
plausible hamiltonian for the interaction of physical qubits. 
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In the following we mention some of the works only in this latter category which are 
of relevance to our work in this paper. 

After the work of Nielson _21_ , it has been shown that two spins interacting by the Ising 
interaction in the z direction, when placed in a magnetic field of arbitrary direction, 
acquire maximum entanglement when the magnetic field is perpendicular to the z 
direction |?7] . 

The effect of anisotropy (in the spin couplings in the x, y and z directions) has also 
been studied in a number of works for different models j281 1291 1^ I31j . The effect of 
inhomogeneous magnetic fields has been studied in |22], but only on an XY system. 
Such a system already shows entanglement when placed in a uniform magnetic field. 

1.3 Results 

In this paper we have studied an isotropic two qubit system in a inhomogeneous 
magnetic field, described by the Hamiltonian 

H = J^^.^2-r{B-rh)a^,-r{B-h)a2., (4) 

where J is the isotropic coupling between the spins, -B > 0, and the magnetic fields on 
the two spins have been so parameterized that h controls the degree of inhomogeneity. 

Let us first review the situation for the homogeneous magnetic field. 

For the ferromagnetic (J < 0) system, there is no thermal entanglement at any 
temperature, but for the anti-ferromagnetic (J > 0) case, thermal entanglement de- 
velops when the temperature drops below the threshold value kTc := We want 
to see how the presence of inhomogeneity modifies this situation. We will show that 
inhomogeneity has the following effects: 

1- In the ferromagnetic system it generally produces entanglement, dependent on 
the value of the magnetic field and the temperature. There is a threshold temper- 
ature above which no entanglement is possible. This temperature has in fact been 
zero in the uniform case which has been shifted to finite values by the inhomogeneity. 
Specially at temperatures near zero and in zero magnetic field, the effect of inhomo- 
geneity is very significant. Under this condition a very small inhomogeneity produces 
maximal entanglement as shown in figures (3 and 4). 

2- In contrast to the ferromagnetic case, the effects in the anti-ferromagnetic sys- 
tem are small. Inhomogeneity in this case slightly raises the threshold temperature, 
and lowers the value of entanglement as shown in figures (5 and 6). 

The structure of this paper is as follows: After presenting the essentials of ther- 
mal entanglement in the next section, in section we study the spectrum of the 
Hamiltonian and characterize the entanglement of the ground state in various regions 
of the parameter space. In section |1] we analyze the thermal entanglement of the 
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system. Throughout the paper we normahze the couphng between spins to J = 1 for 
the anti- ferromagnetic case and to J = —1 for the ferromagnetic case and study the 
results for the two cases separately. 



2 Preliminaries on thermal entanglement 

A spin system with Hamiltonian H kept at temperature T is characterized by a density 
matrix p := ^e~^^ , where P = k is the Boltzman constant and Z := tre~^^ is 
the partition function. 

The entanglement of this density matrix, called the thermal entanglement of the spin 
system can be calculated exactly with the help of Wootters formula 20 . Explicitly 
it is given by the following formula 



E{p) = 2 ^ ^ log2 2 ' 

where 

C = max{0, Ai — A2 — A3 — A4}, (6) 

and A's are the positive square roots of the eigenvalues of the matrix pp in decreasing 
order. The matrix p is defined as 

p = {ay (g)ay)p*{ay (S)ay), (7) 

where * denotes complex conjugation in the computational basis. 
In case that the state is pure p = \^){^\, with 

IV^) ■.= a\+,+)+b\+,-) + c\-,+)+d\-,-) (8) 

the above formula for the concurrence is simplified to 

C{ip) =2\ad-bc\. (9) 

Since E is an increasing function of C, it is usual to take C itself as a measure of 
entanglement whose value ranges from for a disentangled state to 1 for a maximally 
entangled state. In the following sections we apply this formalism to the inhomoge- 
neous system given by the Hamiltonian 



3 Ground state entanglement 

When the magnetic field is uniform, i.e. 6 = 0, the hamiltonian ^ has two sym- 
metries, namely [H, Sz] = [H, S'^] = 0, where Sz and S'^ are the third component of 
spin and the total spin respectively. In a inhomogeneous magnetic field, the sym- 
metry [H, S'^] = no longer holds and thus the triplet and the singlet spins are no 
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longer energy eigenstates separately. A straightforward calculation gives the following 
eigenstates: 



\h) = !-,-> 



'"^'^ = ^2{S' + (1 - 0^; 

with corresponding energies 



^ ((5-1 + 01+, -) + (<5 + l-0|-,+)) 



((5 + 1-01+,-) -(5-1 + 01-,+)), (10) 



El = J + 2B 

E2 = J-2B 

E3 = -J{l-20 

ea = -j(i + 2e), (11) 



where ^ := + 5"^ and 5 = j. 



Note that we are working in units so that B and J are dimensionless. It turns out 
that ^ is the suitable parameter for expressing the effects of inhomogeneity. Thus 
hereafter we will mostly use ^ rather than the original parameter 6, in our analysis. 
The value ^ = 1 corresponds to a uniform magnetic field and deviations from this 
value characterize the degree of non-uniformity. 

In the limiting case — > 1, the two states and |02) respectively go to the 
maximally entangled states :^(l+,— ) + | — ,+)) ^^l"*"'"^ ~ |— ,+))• 

3.1 The ferromagnetic case, J = — 1 

The ground state depends on the value of the magnetic field B and the inhomogeneity 
parameter ^. It is readily found that the ground state energy is equal to: 



^2 = -1-25 if i<B + l, 
E^= l-2i if i>B + l 



(12) 



Thus for ^ < S+1, the ground state is the disentangled state |(/)2) and for ^ > B+1, 
the ground state is the entangled state |i;^3). 
The phase diagram of the ground state is shown in figure (1). 

For each value of the magnetic field B, there is a threshold parameter ^-^ := B + 1 
above which the ground state will become entangled. Conversely for each value of 
inhomogeneity ^ there is a value of magnetic field B^ := ,^ — 1 above which the ground 
state will loose its entanglement. 

In the entangled phase the entanglemenet of the ground state is found from @ and 
(TkHi to be 

cm = J. (13) 
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Figure 1: (color online) The ground state of the ferromagnetic C0jSG, BjS 8i function of 
inhomogeneity ^ and the magnetic field B. We work in units where B is dimensionless. 
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Figure 2: (color online) The ground state of the antiferrmagnetic case, as a function of 
inhomogeneity ^ and the magnetic field B. We work in units where B is dimensionless. 



which is solely determined by inhomogeneity. A very interesting point is that when 
-6 = 0, with an infinitesimal value of 6 ~ 1) the system enters the maximally 
entangled phase |i;^3) with etanglment C = | ~ 1. This remarkable feature means 
that the absence of entanglement in ferromagnetic Heisenberg chain is completely 
unstable against very small inhomogeneities. It is also reminiscent of quantum phase 
transitions where a slight change in one of the parameters of the system, changes the 
behavior of the system dramatically. Increasing further the inhomogeneity will move 
the ground state further into the entangled phase but reduces its entanglement due 
to (p)) . 

3.2 The anti-ferromagnetic case, J = 1 

In this case we find that the ground state energy is equal to 



Thus for ^ < B—1, the ground state is the disentangled state |(/)2) and for ^ > B—1, 




(14) 



the ground state is the entangled state \4>a). The phase diagram of the ground state 
is shown in figure (2). 

Again in the entangled phase the entanglement of the ground state is found from 
(jSJ and (Uni) to be 

C{ct>A) = ). (15) 

which is independent of B. Increasing inhomogeneity again decreases the concurrence 
and hence the entanglement. 



4 Thermal entanglement 



Raising the temperature mixes the ground state with excited states. Depending on 
the sign of J and the value of parameters this may increase or decrease the value 
of entanglement. In some cases the disentangled ground state mixes with entangled 
excited states and in some other cases the entangled ground state mixes with disen- 
tangled excited states. To see what happens exactly we calculate the entanglement 
of the thermal state p = ^e~^^. The symmetry [H,Sz] = constrains the general 
form of p to 



w z 
z w 

V ) 

where C is found from (j^l and I7| to be given ^2 by 



(16) 



C = 2 max (0, 



u 



(17) 



The exact values of the elements of p is obtained by knowing the spectrum of H. 
After a simple calculation from 



(18) 



we obtain 



and 



z 



z = — -e^-^sinh2J/3^, 



(19) 
(20) 
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where Z is the partition function given by 



Z ■= tre-f^^ = 2e"^-^ cosh 2j3B + 26^-^ cosh 2/3 J^. 



(21) 



Thus from (|17() we find that 




(22) 



We consider the ferromagnetic (J = —1) and the anti-ferromagnetic (J = 1) cases 
separately. 

4.1 Ferromagnetic case, J = — 1 

Setting J = —1 in (|22)) we have 



In the uniform case (^ = 1), this equation turns into e^^ = —1 which has no solution. 
Thus in this limit there is no thermal entanglement in the spin system in accordance 
with previous results (HHHHH]). 

However in the inhomogeneous case ^ 1) this equation has nontrivial solutions. 
Figure (7) shows the variation of threshold temperature with ^. 
Figure (3) shows the entanglement as measured by the concurrence for a fixed value 
of inhomogeneity ^ = 1.1 in terms of the temperature and magnetic field. Below the 
threshold temperature (about 0.25 for this value of ^), thermal entanglement develops 
and is maximized for zero magnetic field B. The value of this maximum entanglement 
occurs of course at T = 0, where its value is equal to |, equal to 0.9 in this case. 

Figure (4) shows the value of entanglement in terms of the temperature and the 
inhomogeneity for zero magnetic field. It is seen that at any temperature there is 
a parameter above which thermal entanglement will develop in the system. The 
value of ^0 is obtained from H24|) and increases with increasing the temperature. At 
very low temperatures is very close to 1 which shows that a small degree of inho- 
mogeneity will develop maximal entanglement in the system. 

4.2 Anti-Ferromagnetic case 
Setting J = 1 in (|22)) we obtain 




(23) 



The threshold temperature is obtained from the equation 



e-'^^ sinh2/3e = C- 



(24) 




(25) 
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Figure 3: (color online) Concurrence versus temperature and magnetic field for ^ = 1.1 
in the ferromagnetic system. 
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Figure 4: (color online) Concurrence versus temperature and inhomogeneity in zero 
magnetic field in the ferromagnetic system. 
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The threshold temperature is obtained from the equation 

e2^sinh2/3^ = ^. (26) 

In the uniform case = 1), this equation turns into e^^ = 3 which gives the threshold 
temperature kT^. = . 

In the inhomogeneous case > 1) this equation can be solved numerically, the re- 
sults is shown in figure (7). It is seen that inhomogeneity only slightly increases the 
threshold temperature in contrast to the ferromagnetic case where it had appreciable 
effect. 

Figure (5) shows the entanglement as measured by the concurrence for a fixed value 
of inhomogeneity ^ = 1.1 in terms of the temperature and the magnetic field and 
figure (6) shows the value of entanglement in terms of the temperature and the inho- 
mogeneity for zero magnetic field. Comparing these figures with figure (3) and with 
the corresponding figure of (-23:) we see that in the anti- ferromagnetic case, inho- 
mogeneity has a small effect on the threshold temperature and magnetic field and 
only decreases the value of entanglement once it is developed. Its value is weakened 
by raising the temperature and near the threshold temperature it has a vanishingly 
small eff'ect. It is seen that for any fixed temperature inhomogeneity always decreases 
entanglement, in contrast to the ferromagnetic case. 

5 Discussion 

We have studied the effect of a inhomogeneous magnetic field on the ground state 
entanglement and thermal entanglement of a two spin system. We have shown that the 
effect of inhomogeneity is most pronounced on ferromagnetic spins, i.e. spins coupled 
by ferromagnetic interactions. At zero temperature an infinitesimal magnetic field 
applied to the two spins in opposite directions maximally entangles the two spins. 
It is as if we twist the two spins into an entangled state. This effect also exists at 
higher temperatures but to much less degree. When the coupling of the spins is 
anti- ferromagnetic inhomogeneity can only have a weakening effect on entanglement. 
Although we have derived our results by studying a two spin systems, these results 
may also hold true more or less on spin chains. A parameter like ^ := \J\ + (6^), where 
(6^) is the average of inhomogeneity on all sites, i.e. (5^) := Y^f=i{Bi — {B))"^ , may 
characterize the infiuence of inhomogeneity on the entanglement of a spin chain. 
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